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NOTICES 


When Government drawings, specifications, or other data are used 
for any purpose other than in connection witha definitely related Govern- 
ment procurement operation, the United States Government thereby in- 
cursnoresponsibility nor any obligation whatsoever; and the fact that 
the Government may have formulated, furnished, or in any way supplied 
the said drawings, specifications, or other data, is not to be regarded 
by implication or otherwise as in any manner licensing the holder or 
any other person or corporation,or conveying any rights or permission 
to manufacture, use, or sell any patented invention that may in any way 
be related thereto. 


The information furnished herewith is made available for study 
upon the understanding that the Government’s proprietary interests in 
and relating thereto shall not be impaired. It is desired that the Judge 
Advocate (WCJ), Wright Air Development Center, Wright-Patterson 
Air Force Base, Ohio, be promptly notified of any apparent conflict be- 
tween the Government’s proprietary interests and those of others. 
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FOREWORD 


Computations presented in this report were accomplished by 
Dr. Clarence Ross and Mr. Carl S. Fluke of Computation Branch, 
Aeronautical Research Laboratory, Wright Air Development Center. 
Work was requested by Air Force Cambridge Research Center in 
relation to RDO No. 114-26E, ‘‘Analysis of Geophysical Intelli- 
gence.’’ Dr. W. Pfister of the Geophysics Research Division, 
AF Cambridge Research Center, Project Scientist, has supervised 
and approved solutions presented. 


Equations and measurements used in the computations were 
furnished by Pennsylvania State College under contract AF 19(122)-44. 
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ABSTRACT 


This report gives the results of extensive numerical calculations 
of the index of refraction, the absorption index, and the complex polari- 
zation vector for radio waves of low frequencies being propogated in an 
absorbing ionized medium in the presence of an external magnetic field. 
All of the calculations reported here were performed on OARAC, a new large 
scale digital computer operated for WADC at Wright-Patterson AFB, Ohio. 
The results are presented in tabular form. Any interpolation in these 
tables must necessarily make use of divided differences because of unequal 
argument intervals. | 


PUBLICATION REVIEW 


This report has been reviewed and is approved. 


FOR THE COMMANDER: 


| a7 
Be 
7 Aas B. re 


_f Colonel, USAF 


wv Chief, Aeronautical Research Laboratory 
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~ INTRODUCTION 


This report comprises all of a large scale program proposed by the 
Atmospheric Ionization Laboratory of Pennsylvania State College in June 
1950. Herein are given the results of extensive numerical calculations 
of the index of refraction, the absorption index, and the polarization 
vector for radio waves on a carrier frequency of)from 25 to 500 kilo-. 
cycles per second propagated vertically in the lower ionosphere with the 
magnetic field conditions at State College, Pennsylvania. 


The full proposal consists of two parts. Part A Reaueets calcula- 
tions at ten frequencies, approximately 30 values of collision frequency, 
‘and approximately 35 values of electron density, all parameters to be 
varied separately. Part B contains calculations for.a frequency of 150 - 
kilocycles/second using ten values of collision frequency and eleven 
values of electron density to be repeated for three different values of . 
syromagnetic frequency, and three different values of the dip angle. 
Again all parameters are to be varied separately. This makes a total of 
10,858 sets of parameters. | 


| Since this agency is concerned solely with the computations involved, 
no study of the results or discussions will be given. This will be done 
by the proposing agency. | - 2 


The primary theoretical attack of the proposing Laboratory is direct-— 
ed towards explaining or predicting theoretically the magnitudes and vari- 
ations of the quantities measured there, such as group path, reflection 
coefficient, phase path changes, polarization, and fading rates. Past 
experience with the comparatively simple ray theory has shown that, for 
low frequencies, it is difficult indeed to carry out the integrations 
needed to obtain the derived quantities, even when the assumed ionization 
and collision frequency distributions are rather simple. When the distri- 
butions used are complicated, as would be indicated to be necessary by 
experimental results, the hope of obtaining analytic solutions is consider— 
ably less. 


Moreover, experiment and theory both appear to indicate that ray 
theory is insufficient and that some form of wave theory must be used. In 
its simplest form, this wave theory requires that one obtain the solution of 
a second order differential equation with a variable coefficient. As it is 
well known, this equation cannot, in general, be solved exactly. 
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Therefore, it has been felt that the most fruitful approach for the 
proposing laboratory will involve some form of numerical solution, at least 
insofar as the parameters inserted into the differential equations are con- 
cerned. One such procedure for solving this equation has been developed 
by other members of that group (see Technical Report No. 20, Penn State). 


For this method dealing with the wave equation, as well as for many 
of the more common methods for treating the ionospheric problem, it is 
necessary to know the index refraction, the absorption index, and the 
polarization. The purpose of this study is to supply numerical values 
of these quantities for a wide range of ionospheric conditions for an 
operating frequency of from 25 to 500 kilocycles per second. A complete 
description of these parameters appears later. They were chosen to cover 
a range of values which should be more than sufficient to fit the actual 
physical conditions. All of these calculations were made for magnetic 
field conditions appropriate to vertical incidence at State College, Penn— 


sylvania for approximate level of the E-layer (except for 660 sets ‘in 
Part B). | 


For 1230 of these sets of parameters, the calculations were carried 
out during the actual acceptance test of OARAC by the Air Force. The re— 
mainder of the 10,858 sets were completed after OARAC was installed per— 
menently at this agency. OARAC is a large scale, automatic sequencing, 
digital calculator built for this agency by General Electric at Electronics 
Park, Syracuse, New York. This is the first major problem of comparable’ 
$ize to be completed on it. 
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EQUATIONS 


These calculations were based on the well known Appelton—Hartree 
formula. (For references, derivation, and discussion see: 5.K. Mitra, 
"The Upper Atmosphere", Calcutta, 1948). : 


The two equations involved are: 


(1) MY = (peje)? #4 + < 


2 
oy | 2 
. — eerie eae . eee, Sy aes + 


where 


= index of refraction 

= absorption index 

= r+ js = polarization 

= electron density 

= external magnetic field 

= angle between direction of propagation and external iaanerie field 
= 2uf = angular frequency of wave 

= collisional frequency 

- - F/P? 

= Pr/P% 

= PP, /Po 

= H./me 

P, cos. 0, 

= P. sin @. 

b eunconents of the external magnetic field in and perpendicular 
to the plane containing the direction of propagation and the extern- 
al magnetic field. | 
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il 
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e, m = charge and mass of an electron. 


For the actual calculations we used a procedure given by V. A. Bailey, 
Phil. Mag., 26, 425-458, Oct. 1988, which was obtained by applying complex 
transformations to Equations (1) and (2). For a given P, Pa» 6, v, and N 
we evaluate the following quantities in succession: e 


(83) PE =3.19x10°x N- 

(4) o =2 cos 6, /sin® g 

(5) x = ov/Py, y = oP/Py(1-P6/P*) 
(6) pe =xt sy? | 

(7) & = x/p°, n=-y/pe* 

(82) A= -n'°-1, B= 2 En 


(9) Ce |Va2 + BR 
(10) -F = 2 Iiccea)| 


(41) r= &(1+ F), s = n(1 + $) 

(182) X=- y(% te cos @,) -1, Y = y (4 + r cos 0.) 
(18) L® = +? 

(14) a= Ve BRE, b= 1+ X/L 

(15) p= W1 (a+b)! z= WA (a-p) | 


The required results for the bles of .refraction and the absorption 
-, index are given in Equation (15). The quantities r ands ‘of Equation (11) 
. are the real and imaginary parts, respectively, of the polarization, R, de- 
- . fined in Equation (2). It is sufficient to tabulate R in terms of these 
two quantities. 


When the ratio of B?/A? becomes small, the quantity F in Equation (10) 


cannot be computed directly. It is necessary to expand this quantity in a 
series in yee Then the following equation was used in place of Equation (10): 


(16) F = A ef) wnen A <0 and [dhl < 0.001 
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2 4 | 
and F | [1 . p* , ] when A > 0 and | dtl < 0.02. 


This change was made automatically by the calculating machine. 


As can be readily seen from Equation (11), r ands as well as all 
of the quantities in the following equations, are double valued. One 
‘value is given by using the positive sign for both r'and's, while the 
other value is given by using: the negative sign for both r ands. One 
of these two values corresponds to the "ordinary ray", while the other 
corresponds to the "extraordinary ray". 


‘ fj 
The signs to be used for the components are given in Fig 1, below; 
‘the signs refer to those used in Equation (11). 


a A a 
Eurortiey [> | >| > > 


| : | p 
Where v_ is the critical collision frequency defined by vo = 


and N, is the critical ionization for reflection for ace Qsebeorpbton 
and no external magnetic field, and is given by N, = 3.10 x 10° 


PARAMETERS 


The parameter choices are not uniform, but rather are based on 
anticipated behavior, both of the calculations themselves and of the 
signal in the atmosphere. The values of N are grouped around the points 
N, and Ny where 


N. is the value of N for which Pg = P* (as above) 
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, P 
N, is the value for which Pg = P*(1 + -H). 


In the region of these two points the curves are likely to change rapidly 
when v is small. 


since none of the parameter intervals are equal any interpolations 
in or between the tables must use some form of divided differences. 


The actual parameters used are given with the results, but the 
values of N. and v, are needed in Fig 1 and are given below, as well as 


other important parameters. 


f N. N, Vv. Py * 6. 
Part A 25,000 7.73483 470.900 5.348x10®  9.408x10° 19°8' 
55,000 387.4864 1056.39 " " i 
70,000 60.6408 1857.50 " " " 
80, 000 79.2044 1581.08 " iT " 
100,000 123.757 1976.41 u ; . 
125,000 198.370 2509.20 n t i" 
150,000 278.454 3057.44 " " ". 
200,000 495.029 4200.38 ' , 1 
300,000 1113.8 6671.79 " iN 1" 
500, 000 3093.92 12357 .2 " " T) 
Part B 150,000 278.454 2987.49 3.12077x10® 9.0x10° 15° 
Ty " i" 5. 11715x10° " 19°8' 
tt ie 1 8. 86822x105 " 95° 
" : 3057.44  3.26156x10° 9.406x10° 15° 
: : j 5.348x10° " 19°83" 
| it Uh 9. 2BR27x 1 0 i" 960 
" ‘ 8232.98  3.46758x10° 107 15° 
Li tt i" 5, 68573x10° rT) 19°8' 
it 9 WH 9. 85358x105 " 25 re) 
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CALCULATIONS 


Due to the large range of the parameter values, and the intermedi- 
ate results as well, the entire computations were carried out.using 
multiple accuracy. OARAC is a ten digit fixed decimal computer, but 
certain operations permit the use of twenty or more digits by special 
coding. It was found that twenty digits would give the final accuracy 
desired in most of the computations. However more digits were needed 
in certain places. 


Four subroutines were used to increase the efficiency of coding. 
These were xy, x”, “x and for twenty digits of x and y. These re—- 


| X 
Quired approximately 60% of the two minutes needed to compute each 
case. The total computing time is therefore approximately 345 hours. 


In outscribing the results from magnetic tape, ten digits were 
deemed necessary and sufficient. These were chosen to be three to the 
left and seven to the right of the decimal point. In order to be able 
to print as many significant figures as were accurate, all numbers were 
printed to full ten digits. Occasionally eleven digits were needed to 
tabulate the r,x10* values. Obviously then, all the significant digits 
printed is no indication of the accuracy. 


Hach step in’ computing was checked by duplicating the computation 
by another means. These results were automatically checked and the 
computation was either continued or repeated depending on the outcome 
of each check. Approximately 180 of these checks were coded into the 
problem. The final results are therefore as thoroughly correct as 
possible. The accuracy of the results should be trustworthy to at 
least five significant figures as was desired. Many are more accurate, 
however, depending on the range of the parameters. 


There are 299 sets of collision frequency in Part A and these 
are presented in the Appendix, with one set per page. In part B, there 
are 90 sets of collision frequency and these are rresented following 
. Part A, with two sets per page. 


The values of P,, f, », and @ were tabulated automatically at 
top of each page and the electron density is the running independent 
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variable. Certain values had to be shifted when placed at the decimal 
point for printing and the appropriate power of ten is indicated. The 
"+" and "—" subscripts indicate the use of the sign in Equation (11). 
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